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Context. The hierarchical model of galaxy formation, despite its many successes, still overpredicts the baryons fraction locked in galaxies as 
a condensed phase. The efficiency of supernovae feedback, proposed a long time ago as a possible solution for this so-called "overcooling" 
problem, is still under debate, mainly because modelling supernovae explosions within a turbulent interstellar medium, while capturing realistic 
large scale flows around the galaxy is a very demanding task. 

Aims. Our goal is to study the effect of supernovae feedback on a disk galaxy, taking into account the impact of infalling gas on both the star 
formation history and the corresponding outflow structure, the apparition of a supernovae-driven wind being highly sensitive to the halo mass, 
the galaxy spin and the star formation efficiency. 

Methods. We model our galaxies as cooling and collapsing NFW spheres. The dark matter component is modelled as a static external potential, 
while the baryon component is described by the Euler equations using the AMR code RAMSES. Metal-dependent cooling and supernovae- 
heating are also implemented using state-of-the-art recipes coming from cosmological simulations. We allow for 3 parameters to vary: the halo 
circular velocity, the spin parameter and the star formation efficiency. 

Results. We found that the ram pressure of infalling material is the key factor limiting the apparition of galactic winds. We obtain a very low 
feedback efficiency, with supernovae to wind energy conversion factor around one percent, so that only low circular velocity galaxies give rise 
to strong winds. For massive galaxies, we obtain a galatic fountain, for which we discuss the observational properties. 

Conclusions. We conclude that for quiescent isolated galaxies, galactic winds appear only in very low mass systems. Although that can quite 
efficiently enrich the IGM with metals, they don't carry away enough cold material to solve the overcooling problem. 
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1. Introduction 

The hierarchical model of structure formation, for which 
massive galaxies grow by mergers of smaller satellites or by 
filamentary accretion of gas, is now a well-established theory 
(the so-called CDM paradigm), that compares favorabily with 
observations. In this standard picture, dark matter is the main 
component that drives structure formation via gravitational 
instability from primordial density perturbations. On the other 
hand, numerical simulations including gas dynamics with 
radiative cooling and some recipes of star formation tend to 
produce galactic disks that are too small, to dense, with too 
low angular momentum and too many condensed baryons, 
when compared to o bservations. These v a rious problems, 
called "overcooling" (jPekel & Silk 1986t Blanchard et al 



1998 



2003) 



1992[ Navarro & White! 1 19931: lYepes et al.l ll997b iGnedin 



1999; iCole et al.ll2000t iKay et al.H2002t iSpringel & H ernciuist 



or "angular momen tu m" (Navarro & Bena 1991 
Sommer-Larsen et al. 19991: Steinmetz & Navarro! 1 1999 : 



Bullock et al.ll200lt iMaller & Dekell 120021; lAbadi et al J 12003 



Sommer-Larsen et al. 2003 ; Read & Gilmore 120051) problem, 
may all have the same physical origin. Nevertheless, they are 
considered as weak points in the CDM theory. 

Supernovae-driven winds are a key ingredient of current 
galaxy formation models, in order to suppress the formation of 
low-mass galaxies and maybe to solve the "overcooling" prob- 
lem. They can potentially suppress the fast consumption of gas 
by transferring cold and dense gas into the hot and diffuse sur- 
rounding medium of a galaxy. Thus, a fraction of gas would 
not collapse into cold dense molecular clouds and would not 
transform into stars. That is the most favored explanation for 
the disruption of star formation and the lack of baryons avail- 
able i n stars dSpringel & Hernquistl 2003 : Rasera & Teyssier 
20061; IStinson et alj|200fj|) . Feedback is also invoke d to repro 



iHultman & Pharasvnl 1 19991 ISomerville & Primack duce t he large disc galaxies observ ed at our epoch dWeil et al 
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1998 : Thacker & Couchmanl 200 1]). th e morphology of galac 
tic discs Jsommer-Larsen et al I T20031 lOkamoto et al.l l20o3) . 
to stabilize the gaseou s disc against the Toomre instability 
( Ro bertson et al.l 120041) or to fit the te mperature of the in ter- 
galactic medium (IGM) at redshift ~ 3 dCen & Brvanll200lb . 
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The proper modelling of galactic winds is quite dif- 
ficult to handle, especially in cosmological simulations, 
because of the very large scale separation between the 
galactic outflow in the IGM and the supernova explosions 
within the interstellar medium (ISM). But it represents 
an unavoidable task in order to solve the weakest points 



els dWhite & Frenk 199 it Somervil e & Primackl 19991 


Kauffmann et al. 


19991 Cole et all 200C ; Hattonetal. 


20031 Monaco! 


2004c Bertone et al. 2005) and in numer- 


ical simulation ( Cen & Ostriket 1992c Navarro & White 


19931 Mihos & Hernauistl 1994c Katz et al. 


1996; 


Thacker & CouchmanI 2000c Scannacieco et al. 


2001; 


Sprineel & Hernquistl2003c Rasera & Tevssierl2006l) of £ 


;alaxy 



formation. Observational evidence for galactic outflows have 



already been pointed out by several authors ([Bland & Tully 
1988flHeckman et alJ[l990l:lEUison et al]l200d:lHeckman et al 



2000t iPettimetal l 1200 lL 12002c lAdelberger et alJ 12003c 
Bouche et alJ 20061) . They are usually associated to massive 
starbursts, for which very strong outflows are reported: for one 
solar mass of stars formed in the galaxy, between 1 and 5 solar 
mass of gas ar e ejected in the winds of starburst dwarf galaxies 
dMartinlll999h . 

Despite the difficulties of modelling supernovae explosions 
within a turbulent, multiphase and magnetized ISM, under- 
standing the physics of the re sulting large scale outflows is also 
a challenge. As explored by iFuiita et al.1 ( 12004b in the context 
of an isolated, pre-formed galactic disc, the ram-pressure of in- 
falling material might be the main limiting factor for galactic 
winds to exist. The gravitational potential is likely to play a 
minor role in the development of galactic winds, especially in 
low-mass galaxies where the escape velocity is fa r lower than 



the ve locity of the wind. This fact, outlined first by lFujita et al 
(120041) . is in contradiction with most semi-analytical studies 



so far, for which only the escalpe velocity of the parent halo 
determines if a wind can develop or not. In this paper, using a 
self-consistent modelling of gas cooling and accretion from the 
hot surrounding medium into the disc, we will demonstrate that 
it is indeed the ram pressure of the infalling gas that matters for 
galactic outflows. 

To illustrate further the difficulties in modelling super- 
novae feedback, we note that the mass ejection rates mea- 
sured in galactic winds obtained in numerical simulations are 
widely discrepant. Some authors report very efficient winds, 
even in the quite massive galaxies ( Springel & Hernquistl 
2003c iTasker & Brvan] 120061) . while other authors obtained 
very low ejection effici ency, mainly in dwarf galaxies 



dMac Low & Ferrara Il999h . We believe that these differences 
arise partly from the different numerical technics used, and 
partly from the different boundary conditions imposed around 
the galaxy. In this paper, we report the first self-consistent 
adapti ve mesh refinement (AMR) simulation of a collapsing 
NFW ( Navarro et al. 1996b spherical cloud forming a galactic 
disk from the inside out. In the recent literature, this same prob- 
lem was ad ressed using smooth particle h ydrodynamics (SPH) 
technics by ISpringel & Hernquistl (120031) . Grid-based simula- 



tions were a lso performed, bu t using isolated, pre-fo rmed ro- 
tating discs dFuiita et al.ll2004t ITasker & Brvanll2006l) . One of 



the key outcome of such numerical simulation is the fraction of 
the injected supernovae energy which is transferred up to the 
large scale galactic wind and is truely ava ilable for feedback 
dFuiita et al.ll2004clScannapieco et al.ll2006l) . 

We are indeed facing a major difficulty when introduc- 
ing feedback in numerical simulations that include radiative 
losses: the characteristic time scale of radiative cooling in 
dense regions where star formation occurs is very fast, usu- 
ally faster than the simulation time step, so that one relies on 
implicit time integration schemes. As a consequence, if one 
handles supernovae heating by dumping thermal energy in the 
same dense regions, most (if not all) of the energy is radiated 
away before any large scale flow develops. This spurious ef- 
fect has bee n identified for a long time, both in cosmological 



simul ations dKata 1 1 992l iMihos & Hernquist 11994 Katz et al 
1996[) and in interstellar me dium simulations dWada et al 
2000[|joung & Mac Low|2"006l). We foll ow the recipe proposed 
by lMori et al .1 dl 9971) and ldnedinld 19981) by adding around each 
supernova event density, velocity and energy profiles based on 
the Sedov blast wave solution, in order to get the correct ki- 
netic to thermal energy fraction in the flow. In this way, a 
large fraction of kinetic energy is injected into the disc, with- 
out being radiated away on the spot. This will translate into 
a turbulent gaseous disc, and in some cases into a large scale 
outflow. To account for the small scale effect of supernovae 



feedback, we follow the multiphase approach oflYepes etal 



19971 lAscasibar et all l2002l ISpringel & Hernquistl l2003i and 



Marri & W hite 2003, for which the supernovae thermal energy 



is collecticely acccounted for as a new equation of state for 
the dense ISM, depending on the gas local density (polytropic 
equation of state). In this way, we can reproduce a stiffer equa- 
tion o f state for the dense gas, leading to t hicker, more stable 



discs dToomrdl 19641 iRobertson et al.l l2004). Another problem 



arises because the scales where radiative losses of SNe rem- 
anants occur are largely unresolved (< 1 pc). It is therefore 
likely that th e energy transfer of SN remnants t o the ISM is 



much lower dLarsonlll974t iThornton et alJI 1 9981) than the en 



ergy trans fered from supperbubbles initiated by multiple local- 
ized S Ne dDekel & Silldll986t iMelioli & de Gouveia Dal Pinol 



2004). In reality, modeling properly the energy deposition re- 
quires to resolve the coherence length of the SNe explosions. 
Thus, we have to bear in mind that our supperbubble model 
entails many assumptions. 

Within this framework, we would like to answer the follow- 
ing questions: What is the conversion efficiency between the 
smale scale supernovae luminosity and the large scale outflow 
luminosity ? What are the conditions for a galactic wind to de- 
velop and escape from the parent halo potential well ? What is 
the mass ejection rate of such a wind ? What is the metallicity 
of the wind and other associated observational signatures ? 

Although the cooling NFW halo approach that we present 
in this paper represents an important improvement over previ- 
ous isolated disc grid-based simulations, we are still a long way 
from realistic cosmological simulations. Such idealized simu- 
lations are nevertheless useful to understand the physical pro- 
cesses responsible for the apparition of a galactic wind subject 
to boundary conditions which are relevant for cosmology. They 
also give us the possibility to derive a simple analytical model 
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outer region of the halo are falling into the inner disc. Both 
density profiles are therefore given by 



o.oi 



o.io 



Fig. 1. The epoch of wind formation tw in Gyr as a function of 
the hydrodynamical efficiency X- We have assumed three dif- 
ferent average star formation time < U >= 1 Gyr (thin solid 
line), < U >= 2 Gyr (dotted line) and < U >= 4 Gyr (dashed 
line) for the 1O 1O M halo, and the same < t* >= 3 Gyr for 
the 10 11 M (solid line) and the 10 12 M o (thick solid line) ha- 
los. Open circles are for Sc (% = 3 Gyr) and Sd (tw = 4 Gyr), 
circles are for Sa (tw = 1 Gyr) and Sb (t w = 2 Gyr). 



that highlights the importance of infall in the process of wind 
formation. Our paper is organized as follows. In section [2] we 
present our analytical model and introduce the relevant quan- 
tities for constraining the galactic wind formation epoch. In 
section [3] we present our simulation settings, with the physics 
involved in our different disc simulations, as well as some nu- 
merical details. In section[4] we present our main result, namely 
the formation of either a galactic outflow or a galactic fountain. 
We compare the actual wind formation epoch with our analyti- 
cal predicition. We also describe in details the flow structure in 
each case, and discuss its observational signature. 



2. Galactic winds from cooling NFW halos 

In this section, we will describe the simplified problem we 
would like to study using high-resolution numerical sim- 
ulations, as well as a simple analytical model that will 
guide us in inte rpreting the results. We model an iso- 



P = 



P.v 



r s (\ +r/r s ) 2 ' 
so that the total integrated mass is 

M(< r) = AiTp s r] ln(l + r/r s ) 



1 + r/r s 



(1) 



(2) 



p s is the characteristic density of the halo and r s its characteris- 
tic radius. We define here the Virial radius as the radius where 
the average density is equal to 200 times the critical density, 
assuming Hq = 70 km/s/Mpc. We consider a constant gas frac- 
tion of 15% throughout the halo. The temperature profile is ad- 
justed so that a strict hydrostatic equilibrium of the gas sphere 
is maintained in the non-cooling, non-rotating case. The halo 
is truncated at 2 Virial radii: outside the halo, we imposed a 
constant gas density equal to a small fraction (10 4 ) of the den- 
sity at this radius. Note that this truncation radius is a key pa- 
rameter in this study. The size of the box is taken equal to 6 
Virial radius. The concentration parameter c = r V j r /r s of the 
NFW profile is taken to be 10, independant of the halo mass. 
Initially, our halo is slowly rotating, with an angular momen- 
tum profile corresponding to the average prof ile found in 3D 
cosmological simulations (IBullock et al.ll200ll) . for which the 
specific angular momentum is 



M(< r) 



(3) 



lated NFW halo ( Navarro et al. ( 1996b ). with initially simi- discussed in 



with a spin parameter A = J\E\ 1 ^ 2 /(GM 5 ! 2 ), for which we con- 
sider only two representative values A =0.04 and 0.1. Starting 
from this initial equilibrium configuration, the gas is allowed 
to radiate its thermal energy using standard, metal-dependant 
cooling processes. Note that initially we assume a constant 
metallicity equal to 10~ 3 solar throughout the halo. In agree- 
ment with the standard picture of galaxy formation, a centrifu- 
gally supported disc quickly forms in the halo centre, gradually 
growing from the inside out. Its size depends strongly on the 
chosen spin parameter. When enough dense material is formed 
within the disc, star formation is allowed to proceed in a qui- 
escent way, using a star formation recip e based on a sta ndard 
Schmidt law (see Eq. [TOjl. According to lSalpeteri dl955l) IMF, 
for each solar mass of stars formed, 77^^ — 0.1 solar masses are 
recycled into supernovae ejecta, and drive strong shock waves 
into the surrounding ISM material. 

Our goal is to implement these two rather standard phys- 
ical ingredients in the RAMSES code, in order to study the 
formation of galactic winds. This will be described in the next 
section. In the present section, we would like to develop a sim- 
ple analytical model that allows us to capture the basic mecan- 
ism driving galactic winds. Our picture contrasts with the one 



lar distribution for th e gas and the dark matter components 
( Ascasib ar et al. 2003 ). W e follo w the initial settings proposed 
by ISpringel & Hernquistl (12003b . so our work can be directly 
compared to them. Imposing indentical profiles for the gas and 
for the dark matter is especially important in the outer part, in 
order to have realistic accretion rates at late times, when the 



Fujita et al.l (120041) . who assumed a pre-formed 
gaseous galactic disc, with no self-consistent star formation. A 
blast wave was initiated at the very center of the disc, with an 
energy budget consistent with massive star bursts. Depending 
on the blast wave and disc parameters, two possible scenarii 
emerge: "blow out" for which the disc is completely destroyed, 
and "blow away", corresponding to a galactic wind escaping 
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the galactic disc. These authors have demonstrated that the 
main limiting factor for a wind to appear in the "blow away" 
case is the amount of infalling gas, coming from the surround- 
ing halo. In the present study, we will address this question 
using a different approach, based on quiescent star formation 
and infall from a cooling NFW halo. 

Let us consider that the gas is cooling instantaneously, so 
that pressure support is lost from the very beginning. Each 
spherical shell is therefore free-falling , gradually accelerating 
towards the center. We compute the accretion rate of halo ma- 
terial onto the disc using a standard semi-analytical recipe (e.g. 
Hattonetal 1l2003h . namely by taking the time derivative of the 
integrated accreted mass 



c (< f) = M(< r ff (i)) 



(4) 



where r//(f) is the radius for which the free-fall time f//(r) is 
equal to the current time t. The gaseous disc mass therefore 
grows by accretion and decays because of star formation, lead- 
ing to a simple differential equation for its time evolution 



M g = M acc - 



Mo 



(5) 



Here, < f* > is the (mass-weighted) average star formation rate 
within the disc. It depends on the exact disc geometry (size and 
thickness) and on the Schmidt law parameters. As a first order 
approximation, one can show that 



(6) 



small discs form stars more efficiently than large ones. If one 
assumes that < t* > remains constant in time, one gets the 
following formal solution 



Jo 



M g (t)= | M acc (u) exp(-^— ^-)dn . 

*N f* 



(7) 



We finally deduce the total supernova luminosity in the disc 
using 



Lsn - VsnCsnM* - 7]sn£SI> 



M„ 



< t >* 



(8) 



where €sn — 50 keV is the specific energy produced by one sin- 
gle supernova (10 51 erg for an average progenitor mass of 10 
M Q ). Part of the supernova energy contributes to sustain a tur- 
bulent multiphase ISM. Only a (small) fraction of this energy is 
converted into actual wind luminosity escaping from the dense 
gaseous disc. The fraction of energy that manage to escape 
from the dense gaseous disc will depend on the disc charac- 
teristic (thickness, size, gas content). In our simulations, these 
properties will be specified by the the spin parameter of the 
halo and by the star formation time scale of the Schmidt law. 
We parametrize this unknown conversion factor x by defining 
the wind luminosity as Lw = xLsn- Our goal is to perform nu- 
merical simulations of star forming discs in order to compute 
the actual conversion efficiency resulting from various, com- 
plex hydrodynamical processes within the disc. We thus call 
this parameter ;f the hydrodynamical efficiency. 

Our goal is now to derive under which conditions a galactic 
wind can form within the halo. We know from previous studies 



(IFujita et al.l 120041) that, in presence of infalling gas, the main 
factor preventing the wind from blowing out is the ram pres- 
sure of the surrounding material. This is a much more stringent 
constraint than any criterion based on the escape velocity. If 
we assume that the free-falling cold gas hits the disc plane with 
terminal velocity, we can compute the accretion luminosity as 



c(t) = \vl{r ff (t))M a , 



and 



Voir) 



Jo 



GM(x) 



dx(9) 



We define the wind beak-out epoch tw as the time for which the 
wind luminosity exceeds the accretion luminosity L acc {tw) = 
Lw(tw)- in Figure Q] the wind break-out epoch is shown for 
various halo masses, as a function of our unknown parameter 
X- One clearly sees that the smaller the hydrodynamical effi- 
ciency, the later the wind will blow out of the disc. More im- 
portantely, for each halo mass, there is a minimum efficiency 
below which no wind can form at all. This minimal value for^ 
is around 0.7% (resp. 3% and 15%) for a 10 10 M halo (resp. 
10 11 M Q and 10 12 M ). We conclude from this very simple toy 
model that it is critical to determine the actual hydrodynamical 
efficiency in quiescent star forming galaxies. Our goal is here 
to perform high resolution numerical simulations of isolated 
galaxies, to measure the epoch when galactic winds appear and 
deduce from tw the corresponding (most likely) hydrodynami- 
cal efficiency. 

3. Numerical methods 

Our simulations were performed with the Adapt ive Mesh 
Refinement (AMR) code RAMSES dTevssierll2002l) . The gas 
evolution is computed using a second-order Godunov scheme 
for the Euler equations, while collisionless star particles' tra- 
jectories are computed using a Particle-Mesh solver. The dark 
matter component is accounted for as a constant background 
gravitational potential. Gas cooling is taken into account as 
a source term in the en ergy equation. The cooling function 
is computing using the Sutherland & Dopital ( 1993b cooling 
model, using a look-up table in the temperature and metallicity 
plane. Metals created in the surrounding of supernova explo- 
sions are advected by the hydrodynamics solver as a passive 
scalar. 

We now describe in more details our 3 main physical ingre- 
dients used for the present study, namely star formation, super- 
nova thermal feedback and supernova kinetic feedback. 



3.1. Star formation recipe 

In each cell, gas is converted into star particles following a 
Schmidt law 



p„ = - — if p > po > P* — otherwise , 



(10) 



where po is an arbitrary density threshold defining what we 
consider here as being "interstellar medium". The star forma- 
tion time scale is proportional to the local free-fall time, 



f* — fo 



-1/2 



(ID 
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The two parameters po and to are poorly known and scale de- 
pendent. One standard approach consists in calibrating these 
numbers to t he observed star f ormation rate in local galaxies, 
the so-called iKennicuttl dl998l) law, which translates roughly 
into po - O.lH.crrT 3 and fo - 1 - lOGyr. In the present 
study, we allow to to take 2 values, namely 3 or 8 Gyr corre- 
sponding to star formation efficienci es of respectively ~ 0.05 
and 0.02, last value is compatible with lKrumholz & Tanl ( 2007 ) 
if we suppose that star formation efficiency remains the same 
at very high (> 10 5 H.cm~ 3 ) and unresolved densities (maxi- 
mum density in our simulations stalls to a few 10 1 H.cirT 3 and 
10 2 H.cm 3 respectively for the 10 10 M Q and 10 11 M due to 
polytropic equation of state). When a gas cell is eligible for 
star formation, N collisionless star particles are spawned using 
a Poisson random process, with probability 

P exp(-Ap) , 



P(N) 



Nl 



where the mean value is 
( pAx 3 \ At 



(12) 



(13) 



The mass of a star particle is taken to be an integer multiple 
of m, = poAx 3 )(n /(l + t]sn + t]w) (see Appendix lAl for details 
on T]w), an d is therefore directly con nected to the code spatial 
resolution dRasera & Tevssierll2006l) . The minimum mass of a 
single star particle is therefore ~ 500 M , but high density re- 
gions can spawn more massive star particles. We also ensure 
that no more than 90% of the gas in a cell is depleted by the 
star formation process. 

3.2. Polytropic equation of state 

The second important ingredient for our galactic disc simu- 
lations is the thermal feedback o f supernovae into the ISM, 
based on the multiphase mo del of lAscasibar et al.l (12002b and 
Springel & Hernquistl (12003b . The idea is to assume that the 
ISM is driven by small scale effects (turbulence, thermal in- 
stability, thermal conduction, molecular cloud formation and 
evaporation) that quickly reach steady-state and lead to a 
quasi-equilibrium thermal state, for which the average temper- 
ature is a function of the mass density alone. In practice, when 
the gas density exceeds the ISM threshold po, the gas tempera- 
ture is forced to be equal to or higher than 
yo-l 

(14) 



where jo is the polytropic inde x, and can be chosen to the ex act 
value of the multiphase model ( Springel & Hernquistll2003 ) or 
to a constant value for sake of simplicity. If yp is chosen be - 
tween 4/3 et 2, as prescribed by ISpringel & Hernquistl (12003b . 
the resulting equation of state gives similar results for the disc 
thickness and maximum density. Our standard choice is here 
7o = 5/3. 

3.3. Blast wave model 



The main problem of the previous thermodynamical approach 
is that feedback is only accounted for at mesoscopic scales. On 



the other hand, it is now an observational fact that very large 
supernovae-driven bubbles appear quite naturally in star form- 
ing galaxies. These superbubbles are built by the coalescence of 
numerous supernovae remnants, who collectively create these 
large-scale features, driving very large, macroscopic flows that 
eventually give rise to galactic winds or galactic fountains. 
Current cosmological simulations do not have the necessary 
resolution to resolve individual supernovae blast wave. On the 
other hand, they do have the resolution to resolve superbubbles. 
Our goal is here to inject kinetic energy in the form of spher- 
ical blast waves of size comparable to galactic superbubbles, 
namely rsN — 100 to 200 pc. This kinetic energy injection can 
be considered as a turbulent forcing term with injection scale 

r SN- 

It is now well-established that supernovae energy can 
not be released only by direct thermal energy injection 
( Navarro & Whitelll993 ). For high density star forming region, 
the gas radiates away all this thermal energy in one time step 
due to very fast atomic cooling. The consequence is that SNe 
have no effect on the dynamics. Various methods have been 
proposed to correctly incorporate the SNII kinetic feedback 
into nu merical simulations. Th e first idea, discussed for ex- 
ample in lGovernato et al. I d2006h is to artificially stop radiative 
cooling in the region where the SNII explosion occurs, for a 
time long enough for the blast wave to develop and expand. 
This trick allows for a large enough fraction of the SNII energy 
to be converted into actual gas kinetic energy. The second ap- 
pro achjsJXMfirecJi^rrject kinetic energy into the surrounding 



gas. ISpringel & Hernquistl (120031) used for that purpose high 



velocity collisionless particles that travel over a rather large 
distance (up to a few kpc) before being re-incoporated in the 
fluid as new SPH particles. 

In this paper, we direclty add to each flow variable (den- 
sity, momentum and total energy) a spherical blast wave solu- 
tion of radius r$N = 150 pc, assuming that each supernovae 
releases half of its energy in this kinetic form (the other half 
is accounted for the polytropic equation of state). As soon 
as this distance is large enough compared to the grid resolu- 
tion, we minimize spurious energy losses. Each time a "star" 
particle is created, we simultaneously remove an equivalent 
gas mass from the cell that is entrained by the SNII blast 
wave. Recalling that m* is the mass locked into long-lived 
stars, the total mass removed from each star-forming cell is 
m»(l + rjs n + r/w), where the entrainem ent parameter rjw is de- 
fined as in Springel & Hernquistl ( 2003 ). We have used rjw — 1 



in order to reproduce roughly a Sedov blast wave solution in 
typical galactic discs for which the density is slightly greater 
to po. Note that in order to implement this numerical scheme, 
we have used "debris" particles to describe the Sedov blast 
wave profile around each SNII explosion. Details are given in 
the Appendix. The most important parameter is r$N, the ra- 
dius of the spherical blast wave. T his sets the injection scal e 
of the turbulent cascade in the disc dJoung & Mac Low 120061) . 
This val ue, comparable to the ty pical size of galactic super- 
bubbles ( McKee & Ostrikerlll977 ). is however larger than indi- 
vidual supernovae remmants, whose size just before the snow 
plow phase varies from 10 to 50 pc, depending on the exact 
early shock wave dynamics, whether adiabatic or evaporative 
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Fig. 2. Cut of the gas temperature in the Oyz plane for the Sd simulation at different epochs. The two top pannels are a zoom 4 
times of the simulation box and the two bottom pannels are the entire simulation box. The colour scale gives the temperature in 
log(K). Note that length scales are not the same in each pannel. 



( Cioffi et al.ll988 ). We have not tried to study the impact of this 
important parameter on our results, in order to avoid clouding 
the purpose of this paper. We rather keep this parameter fixed 
to its fiducial value of r SN - 150 pc dMcKee & Ostriker|[l977 , 
Shull & SiHd[T979h . 



3.4. Simulation parameters 

We use for the fluid solver simple "outflow" boundary con- 
ditions, for which all flow variable are assumed to have zero 
gradient at the boundary. All simulations are initialized with a 
NFW gas and dark matter halo in hydrostatic equilibrium, with 



the same density profile for both fluids. The dark matter com- 
ponent is not explicitly simulated here; only baryons (stars and 
gas) are self-gravitating. We therefore add a fixed analytical 
gravitational potential to the solution of the Poisson equation. 

We have performed our simulations for 2 different halo cir- 
cular velocity, namely V V! > = 35 and 75 km/s, which corre- 
spond to Virial masses of 10 10 M and 10 1 1 M Q at redshift zero. 
We have explored two different star formation efficiency, with 
timescales to — 3 Gyr and to — 8 Gyr. We use a coarse grid of 
128 3 cells, in order to properly describe the large scale dynam- 
ics within the halo, which corresponds, in the RAMSES termi- 
nology, to a miminum level of refinement l m - m = 7. The grid 
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is further refined up to 4 or 5 additionnal levels of refinement 
({max = 1 1 or 12 corresponding to ~ 75 pc, see tableQ]for more 
details), based on a quasi-Lagrangian strategy, for which a cell 
is refined if its mass exceeds a threshold of 2 x 10~ 6 the Virial 
mass. This refinement strategy verifies the Jeans length crite- 
rion inside the galactic disc because of the polytropic compo- 
nent that helps to stabilize against fragmentation. Starting with 
2x 10 6 cells, our simulations usually reach a total of 4 to 5 x 10 6 
cells. We end our simulation after a typical integration time of 
6 Gyr, which corresponds to the free-fall time of the halo outer 
boundary at 2 Virial radius. 



4. Simulation results 



In this section, we present our simulation results for the 2 dif- 
ferent halo mass considered in this paper, namely 10 and 
10 11 M . In the first case, a strong wind generally develops af- 
ter a few Gyr. We will examine in details the wind structure 
and physical properties. For the second case, because of the 
ram pressure of infalling gas from the surrounding halo, we 
see no wind forming: a clear galactic fountain sets in and cre- 
ates a turbulent, hot atmosphere above (and below) the rotating 
disc. We will study the structure of this galactic fountain, and 
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Table 1. Parameters used for the simulations performed in this 
paper. 



compare the various density and velocity profiles obtained in 
both cases (wind versus fountain). 



4.1. 10 1U M Q halo 

In figure [4] we have plotted the star formation rate (SFR) ob- 
tained in the low mass galaxies (Sa, Sb, Sc and Sd, see TableQ] 
for details). The smallest disc with the highest star formation 
efficiency (run Sa) shows, as expected, the highest peak in the 
star formation rate. Using our analytical model, we are in a po- 
sition to roughly recover the same curve using f* 1 Gyr. 

Run Sb and Sc have roughly the same value for the prod- 
uct toAo, and therefore, as predicted by our simple analytical 
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2003[). This also corresponds to the "superwind" geometry 



t (Gyr) 

Fig. 4. Star formation rates of the 10 10 M G halos with A = 0.04 
and f = 3 Gyr (solid line), A = 0.04 and t = 8 Gyr (dotted 
line), A = 0.1 and to = 3 Gyr (dashed line) and A = 0.1 and 
to = 8 Gyr (dash-dotted line). 



model (Eq. |6]|7]), have roughly the same star formation history 
(f* ^ 2 Gyr). The largest disc with the lowest star formation 
efficency (f* ^ 4 Gyr) has the lowest SFR. At late time, after 3 
to 4 Gyr, the 4 curves converge to roughly the same value. This 
corresponds to the late epoch when star formation is mainly de- 
termined by the accretion of gas infalling from the outer halo, 
as predicted by equation[7] 

Figure [2] shows a map of the gas temperature in a cut 
perpendicular to the galactic plane for run Sd (1O 1O M with 
A — 0.1 and to = 8 Gyr) at different epochs. The first pannel, 
at t = 1 .5 Gyr, shows no sign of wind, although it corresponds 
already to the peak of star formation in this galaxy (around 0.3 
M /yr). We see several hot features trying to break out of the 
cold disc, but none of them is able to expand significantly, be- 
cause of the ram pressure of infalling halo material. The second 
panel corresponds to a later epoch when the star formation rate 
has dropped by 50%, down to 0.2 M /yr (see also Figure [3). 
Nevertheless, since gas accretion has also dropped signifi- 
cantly, two hot and large cavities now emerge out of the disc 
and propagate in the halo. At t — 4.5 Gyr the galactic wind has 
finally escaped the halo outer boundaries, and is now in a quasi- 
permanent regime. This 2 large cavities have created some kind 
of tunnel, facilitating the ejection of hot, supernovae-driven 
material from the inner disc, all the way out to the intergalactic 
medium. At t = 6 Gyr, we see a final snapshot of the galactic 
wind, with a typical noozle-like shape (see Fig. |2]and[3]i. The 
hot, metal-rich outflowing gas is surrounded bu a dense shell of 
compressed and cooling halo gas. This dense shell will eventu- 
ally fragment into small clouds that will fall back to the disc. 
We believe that this dense shell corres ponds to the structure 
seen in many observations around M82 (IHeckman et al. 1990: 



described in iTenorio-Tagle & Munoz-Tunonl ( 1 19981) (see also 
Veilleuxetal.112005 and references therein). 



4.2. Wind formation epoch 

In order to describe the outflow more quantitatively, we define 
the net mass flux across a shell of radius r mm and thickness 

— r max 1~ mi n aS 



pu.n 



4nr 2 dr = F++F-, 



(15) 



where n-rjr. This net flux is further splitted into two different 
contribution: the positive flux, F+, corresponding only to out- 
flowing volume elements, namely those that satisfy u.n > 0, 
and the negative flux, F_, corresponding to inflowing volume 
elements, satisfying u.n < 0. We focus here on the positive 
flux, in order to detect the outflow, and to estimate the amount 
of gas expelled as a function of time. The shell thickness was 
set to Ar = 2r s , and we vary the shell radius from 5r s to 20r,, 
in order to define a proper radius where to detect the outflow. 
Figure [5] shows the positive flux measured at 4 different radii 
for our 4 low mass galaxies, as a function of time. The wind 
is detected when the positive flux sharply rise from zero to 
its maximum value. The smaller the shell radius is, the earlier 
the wind is detected. For our smallest shell radius, r ml „ — 5r s , 
the measured mass flux remains quasi constant in time, while 
for larger shell radii, the measured mass flux sharply rises to 
a much larger value, namely 2 to 3 times the inner flux, then 
slowly decays to the correct wind mass flux. This means that 
when the wind breaks out of the halo, a significant fraction (up 
to ~ 2/3) of the wind is filled by the hot gaseous halo and ex- 
pelled out of the halo boundaries. When the permanent regime 
sets in at later time, the mass flux is roughly equal to the inner 
one, independant of the radius of the shell. We therefore con- 
sider only the case r mm - 5r s as the correct proxy for measuring 
the mass outflow rate in the wind. 

We have plotted the positive flux at r„„„ = 5r s for our 
4 galaxies in figure [6] The highest SFR galaxy, Sa, corre- 
sponds to our earliest wind formation epoch, around t w \Gyr. 
According to our analytical model (see figure [TJ, this corre- 
sponds to the hydrodynamical efficency parameter ;f 0.9%. 
Although runs Sb and Sc have identical star formation his- 
tory (see figure HJ, they have different wind formation epochs, 
respectively t w - 2 and 3 Gyr. This demonstrates that for a 
given star formation rate, a more compact disc is more efficient 
in forming galactic winds. According to our analytical recipe, 
run Sb corresponds to^ 0.8%, while run Sc barely reaches 
X — 0.6%. The final galaxy, Sd, features the latest wind to de- 
velopp, with t w 4 Gyr, corresponding also to a slightly lower 
hydrodynamical efficency x - 0.7%. The main conclusion of 
this section is that the conversion efficency of supernovae en- 
ergy deposition rate into wind luminosity is rather low in this 
model, around^ 1%. This efficiency is directly related on 
the choice of supernovae energy injection rate, a quite uncertain 
parameter depending on the scale injection energy (here chosen 
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Fig. 6. Mean flux of mass outflowing calculated between r - 
[5r s ;7r s ] for the 1O 1O M halos with A = 0.04 and to = 3 Gyr 
(solid line), A = 0.04 and to = 8 Gyr (dotted line), A - 0.1 and 
fo = 3 Gyr (dashed line) and A — 0.1 and fo = 8 Gyr (dash- 
dotted line). 



Fig. 7. Mean metallicity of the wind calculated between r — 
[5r s ;lr s ] for the 1O 1O M halos with A = 0.04 and t = 3 Gyr 
(solid line), A = 0.04 and to = 8 Gyr (dotted line), A - OA and 
fo = 3 Gyr (dashed line) and A = 0.1 and to = 8 Gyr (dash- 
dotted line). 



equal to 50%). We can see directly from figure Q] that, accord- 
ing to our numerical experiments, no wind will ever form in 
higher mass halo, as will be demonstrated in the next section. 

4.3. Wind efficiency and metallicity 

Before describing the simulation results for the 10 11 M Q halo, 
we would like to characterize the physical properties of the 
galactic winds we have obtained. It is common to define the 
wind efficency as j] w = F+/M*. In other words, the mass out- 
flow rate is expressed in units of the global star formation rate. 
From figure |4] and figure [5] we can compute direclty tj w : we 
find the maximum efficency at late time, when the permanent 
regime is settled, with values ranging from j] w 0.05 for Sc 
and Sd, up to j] w ^0.1 for Sa and Sb. Here again, it appears 
clearly that more compact discs (low spin parameter Aq) gives 
more efficient winds, while extended discs (high spin parame- 
ter) are less efficient. All these values are ex tremely low co m- 
pared to the wind efficiencies observed by Martini ( 19991) in 
Lyman break galaxies with values greater than 1 . These high- 
redshift galaxies are likely massive starbursts, for which our 
quiescent approach does not apply. On the other hand, we have 
determined that compact star formation sites result in earlier 
and stronger winds. Nuclear starbursts can be natural candi- 
dates to host very efficient winds with rj w > 1, although these 
bright galaxies are probably rather rare objects in the early uni- 
verse. 

As seen from figure|3] the hot gas carried away in the wind 
is highly enriched with metals (Z M , Z ). We define the aver- 
age wind metallicity as the positive flux of metals divided by 



the positive flux of mass Z w = (pZu)+/ F+. Figure Q shows the 
mean metallicity of the wind for our various low mass galaxies. 
We see no clear trend of the wind metallicity as a function of 
disc size or star formation efficiency. We can only observe that 
the metallicity of the wind is of about 0.5 to 1 Z Q . This rather 
high value confirms that the wind comes mainly from the gas 
within the galactic disc that is directly enriched by exploding 
supernovae. 

We now briefly discuss how these winds might explain the 
observed metallicity of the intergalactic medium (IGM). We 
know from our simulations that the typical mass outflow rate is 
around M H , 0.01 M /yr (see figure|6]l, with an average wind 
metallicity of Z„, 1 Z and a typical wind velocity around 
u w — 300 km/s (see figure [T3l . If we assume that the IGM 
is photo-ionized to a temperature around 10 4 -10 5 K, we can 
compute the volume V occupied by the wind when the bubble 
stalls by pressure equilibrium with the IGM as 

«> 2. PigmV 



7-1 



(16) 



assuming the gas is adiabatic. After 6 Gyr, the radius R V 1 ^ 
of the expanding wind lies around 1 Mpc. This corresponds to 
a swept up mass of baryons roughly 5 x 10 9 M . If we assume 
that the metals carried out by the wind perfectly mix with the 
zero-metallicity IGM, we obtain for the metallicity of the IGM 
the following very simple relation 



'IGM 



^1 Z, 

u 



'Q 



10- 



io- 2 z 



(17) 



where c s 10-30 km/s is the sound speed of the IGM. Note 
that this analysis is valid only for the idealized case we have 
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Fig. 8. Star formation rates of the 10 11 M G halos with A = 0.04 
and f = 3 Gyr (solid line), A = 0.04 and f = 8 Gyr (dotted 
line), A = 0.1 and fo = 3 Gyr (dashed line) and A = 0.1 and 
fo = 8 Gyr (dash-dotted line). 



considered here. Observations of the Lya forest clouds show 
that they are metal enriched from Z ^ 10~ 3 up to 10" 2 Z o 
dSongaila & Cowid 1 19961 Ellison et al. 20001) . Cosmological 
simulations also confirm this trend, although the various wind 
models proposed (including our) m ight be difficult to reconcile 
dGnediiJll998l,lThacker etaU2002l) . 



4.4. 10 n M Q halo 

We have plotted in figure [8] the star formation history for the 
4 galaxies hosted by the 10 11 M halo. They are qualitatively 
very similar to the 10 10 M Q case, with La and Ld as the most 
extreme models, and Lb and Lc being quasi undistinguishable. 
Each galaxy has however a peak star formation rate which is 
50% higher than its rescaled 1O IO M counterpart. The circu- 
lar velocity in these new galaxies is higher than the low mass 
ones, while the gas sound speed, because of the polytropic mul- 
tiphase model, remains roughly constant, around 10 km/s, so 
that the disc is now much thinner. This results in higher total 
star formation efficencies. Using our analytical model, we can 
reproduce the 3 star formation histories using for run Sa, Sb 
(and Sc) and Sd the values f» ^ 0.75, 1.5 and 3 Gyr respec- 
tively. 

The other striking differences with low mass galaxies is the 
clear presence of an accretion shock surrounding the disc, and 
the notable absence of galactic winds. We observe instead a 
nice galactic fountain, with plumes of hot gas rising above the 
disc, cooling down and falling back to the disc (see figure[9]l. 

The formation of an accretion shock above a certain halo 
mass is now well understood: t he stability analysis of a r adia- 
tive shock wave, performed by iBirnboim & Dekell (12003b . ex- 



plains why for low mass halo no accretion shock forms. The 
mass threshold above which the shock forms depends on the 
metallicity of the cooling gas. We see from figure [9] that in 
our simulation the shock h eated gas above and below the disc 
is metal poor: in this case, IBirnboim & Deke] (2003) found a 
critical mass for shock appearance slightly below 10 11 M for 
a zero metallicity gas, in good agreement with our numerical 
experiment. The accretion shock converts the infall kinetic en- 
ergy into internal energy, so that our analytical treatment of 
wind formation remains valid, replacing the ram pressure by 
the thermal pressure. The accretion shocks reaches a maximum 
radius of ~ 75 kpc what is ~ 2/3 r t ,, r quite the same virial ra- 
dius fraction than the one found in Birnboim & DekeJ ( 2003 ) 
for the cooling case after ~ 4 - 5 Gyr, and finally contracts 
down to ~ 50 kpc at later time. Although the disc is strongly 
perturbed by supernovae blast waves, no galactic winds can 
form, in agreement with our simple analytical prediction, if we 
assume for the hydrodynamical efficency^f ^ 1%, a value close 
to the one obtained in the 10 10 M case. The atmosphere above 
and below the disc is however highly turbulent and perturbed 
by buyoantly rising plumes. This convective flow is more visi- 
ble in the metallicity map of figure [TOl 

Metals are confined within a rather small distance relative 
to the disc plane, namely ~ 20 kpc. In this case, the thermal 
pressure of the shock-heated halo gas is very efficient in con- 
fining the galactic wind. Figure QT| compares the metallicity 
profiles of the gas between the 1O 1O M halo and the 10 11 M 
halo (run Sd versus run Ld) at different epochs. We define for 
that purpose a mass-weighted average metallicity as a function 
of height and a mass-weighted metallicity as a function of ra- 
dius as 



C"pZ2nrdr f + !'" pZ2ndz 

z(z) = p ; z and z(r) = X ; : 



(18) 



where p is the gas density and Z is the gas metallicity. The 
disc radius r</ is taken equal to r s and the disc scale height 
hd = In the wind case, metals are entrained at very 

large distance above and below the disc in the wind noozle. 
The metallicity is roughly solar and tends to saturate at that 
value. In the fountain case, the metallicity in the halo remains 
rather low, roughly one tenth solar. On the other hand, the gas 
metallicity in the central part of the disc is significantly higher 
than in the wind case (Z gas 1 .2 Z versus Z gas 0.6 Z , see 
figure ITTb. This could be interpreted as the effect of the wind 
removing metals preferentially from the bulge of the galaxy. In 
the outer parts of the disc, however, the gas metallicity is very 
similar in both cases. This could also be interpreted as a higher 
star formation efficiency in the high mass case, especially in the 
central p arts where the dens ity is higher. As explained in great 
details in Dalcantonl ( 20061) . it is quite difficult to disentangle 
the influence of star formation, gas infall and winds on metal 
enrichment in galaxies. We believe that, even in our idealized 
case, the metal content in our galaxies is determined by a subtle 
balance between these 3 processes. 
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Fig. 9. Cut of the gas temperature in the Oyz plane for the Ld simulation at different epochs. The 4 pannels are a zoom 4 times of 
the simulation box. The colour scale gives the temperature in log(K). Note that length scales are not the same in each pannel. 



4.5. Stellar distribution and metallicity 

Figure[T2]shows the surface density profile of stars and the cor- 
responding mass-weighted average metallicity as a function of 
the disc radius. In both cases, we see a clear exponential disc, 
with a bulge-like feature in the central part and a sharp de- 
cline of surface density in the outer part. We find a bulge ra- 
dius of roughly 1 kpc and a disc exponential radius of 2.5 kpc 
for the low mass system and of 4 kpc for the high mass one. 
We note that the outer radius of the low mass case (8 kpc) in 
smaller than the rescaled outer radius of its high-mass coun- 
terpart (10 kpc). This is due to the different disc thickness, so 
that for a given surface density, the threshold for star forma- 



tion is reached more easily by the most massive galaxy. This 
confirms the fact that the polytropic multiphase model has a 
great importance in regulating the gaseous disc properties and 
therefore star formation efficiencies. The break in the stellar 
surface density corresponds roughly to the /it average radius 
where the gas reaches the density threshold of star formation 
{po — 0.1 H.crrT 3 ). Note that the overall amount of stars is very 
similar in both cases, demonstrating that the presence of the 
galactic wind has no impact on the amount of baryons locked 
into long-lived stars. This is consistent with the rather low 
wind efficiencies we have measured (tj w 10%). Our simu- 
lations will not solve the overcooling problem. More impor- 
tantly, the stellar metallicity seems also to be rather insensitive 
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Fig. 10. Cut of the gas density (up) and metallicity (bottom) in the Oyz plane for the Ld simulation at different epochs. The two 
left pannels are a zoom 4 times of the simulation box and the two right pannels are the entire simulation box. The colour scale 
gives the density in log(ctrT 3 ) (up) and the metallicity in Z (bottom). Note that length scales are not the same in each pannel. 



to the "wind versus fountain" scenario. We have Z* 0.4 Z 
in the central part of our 10 11 M Q halo galaxy, while we have 
Z* 0.3 Z in the low mass case. Note that we have computed 
a mass-weighted stellar metallicity so we are more sensitive to 
the overall star formation history rather than to the latest OB 
stars formed. Nevertheless, our wind model has no dramatic 
impact on the simulated stellar population, although there is a 
weak trend for low mass galaxies to have slightly s maller stel- 
lar me tallicities than high mass. Here again, as in iDalcanton 
d2006l) . it is quite difficult to conclude whether this trend is due 
to a wind or to a smaller star formation efficency. 



4.6. Gas kinematics 

The most striking and observable differences between the wind 
and the fountain can be found in the gas kinematics in the halo 
atmosphere, below and above the disc. In order to mimic as 
much as possible the observational signatures of our simulated 
galaxies, as seen by a spectral lines analysis from a distant ob- 
server whose line of sight is aligned with the disc rotation axis, 
we define the average metal velocity as 



Vx(z) = 



C * pZv z 2nrdr 
pZ2nrdr 



and v]_{z) 



L d pZv^lnrdr 
C d pZ2nrdr 



(19) 
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Fig. 11. Mean metallicity of the gas as a function of the projected height (left) and radius (right) for Sd simulation (up) and the 
Ld simulation (bottom) at different times: t — 3 Gyr (solid line), t = 4.5 Gyr (dotted line), t — 6 Gyr (dashed line). 



The velocity dispersion and the sound speed are defined as 
usual by 

to to [''' pZc 2 s 2nrdr 

ai(z) = v 2 ± (z) - vAzf and c](z) = ^ (20) 

J pZlnrdr 

Figure Qj] shows these various profiles as a function of height 
for the wind case (left) and the fountain case (right). In the 
former case, we clearly see a strong symetrical expansion of 
the metals, with a strong acceleration up to 200 km/s in 10 kpc, 
then a plateau afterwards when the wind reaches its terminal 
velocity. In the latter case, on the contrary, we see a converging 
flow, with a velocity of the order of 100 km/s, similar to the 
velocity of infalling halo material. At an altitude of 20 kpc, this 
bulk infall velocity declines and in the same time, the velocity 
dispersion rises: we are entering the galactic fountain for which 



the flow is subsonic and highly turbulent. Note that in the wind 
case, the velocity dispersion also rises up to 100 km/s as we 
approach the disc. This can't be attributed to a turbulent flow, 
since it would have been sueprsonic (the sound speed is around 
50 km/s). It is in fact due to a geometric effect: as soon as the 
wind expansion velocity remains subsonic, below an altitude of 
10 kpc, the section of the wind increases with height. When the 
wind turns supersonic, around z — 10 kpc, the section of the 
wind is forced to remain constant. This process is responsible 
for the typical noozle-like shape of the wind. At low altitude, 
since the wind is expanding quasi isotropic ally, this mimics a 
supersonic velocity dispersion that vanishes at higher altitude 
when the wind flow is quasi parallel to the galaxy rotation axis. 

The 10 10 M Q halo is characterized by a fast supersonic flow, 
expanding quasi-isotropically at low altitude but quasi-parallel 
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Fig. 12. Surface density (left) and mean metallicity (right) of the stars as a function of the projected radius for the Sd simulation 
(up) and the Ld simulation (bottom) at different times: t — 3 Gyr (solid line), t = 4.5 Gyr (dotted line), t — 6 Gyr (dashed line). 



to the rotation axis at higher altitude. The 10 n M Q halo, on 
the contrary, is characterized by a converging pattern quickly 
thermalized close to the disc in a subsonic turbulent flow. 



5. Conclusion 



We applied the method used by iProchaska et al.l (120071) to 
test whether our simulations show velocity dispersion in the 
metal absorption lines of QSOs wider than the velocity dis- 
persion expected from pure gravitational processes (of the or- 
der of the galaxy circular velocity). In this way, we can test 
whether our winds wo uld have been detected. In their paper, 
Prochaska et al.l (12007b use AV90 as the velocity interval encom- 
passing 90% of the mass of a given metal absorber. We simulate 
an observer watching the galaxy along the direction of the wind 
and found that 90% of the mass of metals corresponds to met- 
als with velocity lower than 5 km.s -1 . This is in agreement with 
our findings that wind efficiencies are very low in our model. 



In comparison with ISpringel & Hernquisa (120031) . we have 
very similar initial conditions but a lower wind for r\ w . We 
have therefore more difficulties to solve the 'over-cooling' 
probl em. The explanation has tw o reasons : on the one 
hand ISpringel & Hernquisti (120031) model supernovae explo- 
sions with a phenomenological apporach so that every single 
ejecta in the disc entirely participate to the large scale out- 
flow. In this work, our goal was to model supernovae explo- 
sions with a self-consistent Sedov blast approach where ejec- 
tas are tighlty coupled to the galactic disc (injecting energy in 
the galactic turbulent cascade). On the second hand, we have 
a cooling function depending on metalli city: this increases ra- 
diative losses significantly com pared to ISpringel & Hernquist 
(12003b . Kobavashi et al.l (120061) performed SPH simulations of 
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isolated disks with star formation and self-consistent super- 
novae feedback, metal dependent cooling and the same set of 
initial conditions than ours. They observed very high mass ejec- 
tion rates for the 1O 1O M halo with a total mass ejected ris- 
ing to 80% of the total baryons mass. They noted however a 
strong dependance of their results to numerical resolution and, 
more importantly, to the size of the feedback radius (it can be 
as high as 10 kpc ! ). They even ob s erved a faint wind in the 
10 12 M Q halo case. Tasker & Bryan d2006l) have simulated an 
pre-formed isolated disk with the AMR code ENZO with star 
formation and self-consistent supernovae feedback released in 
pure thermal form over one dynamical time scale. They suc- 
ceed in shuting down star formation but they didn't tak e into 
account infalling gas from the halo. iFuiita et alj ( 120041) con- 
structed an analytical prescription to constrain the mass of the 
halos that can form large-scale outflow with starburst physics. 
They have performed several sets of 2D grid simulations of a 
pre-formed galaxy with an imposed analytical infa ll model and 



found results very similar to ours (see Figure 1 in iFujita et al 



20041) : the key role played by infall and a very low hydrody- 
namical efficiency for quiescent modes of star formation. 

Using a quiescent model of star formation in isolated galax- 
ies, self-consistently simulated from a cooling NFW halo, we 
have studied the conditions for a galactic wind to break-out 
of the ram-pressure exerted by the infalling gas. Our simula- 
tions have shown that a galactic super wind forms in halos of 
10 10 M o and no wind can form in halo of mass greater than 
10 11 M Q , even for our most favorable couple of halo parame- 
ters (A = 0.04, to — 3 Gyr). For large galaxies, a galactic foun- 
tain appears and expells metal rich gas that cannot escape from 
the galaxy. Those results are consistent with our toy model that 
predicts no wind in those galaxies if the hydrodynamical ef- 
ficiency is of a few percents as suggested by our simulations. 
Using this rather simple toy model, we understand this failure 
as due to the ram pressure of infalling material confining the 
outflowing wind. 

Although galactic winds develop in 10 10 M halos, they are 
not sufficient to explain the "overcooling" problem. The mass 
ejection efficiency obtained in our simulations ~ 3 — 12% is 
one or der of magnitude lower than that observed by Martini 
19991) in starburst galaxies. The main conclusion of this pa- 



per is that with a self-consistent treatment of supernovae feed- 
back in idealized simulations of quiescent star-forming galax- 
ies we cannot reproduce the observed mass ejection rates and 
therefore reduce the cold baryon fraction. On the contrary, 
supernovae feedback is an efficient process for metal enrich- 
ment of the IGM. Cosmological simulations with supernovae 
feedback could test this scenario in a more quantitative way. 
Using a more realistic cosmological setting may result into 
non-spherical accretion flows, and therefore to a less stringent 
criterion for a wind to break-out. 

Considering that small mass galaxies preferentially forms 
along filaments, non-spherical accretion should allow the wind 
to break-out more easily and probabily to be more efficient. 
For more massive galaxies, however, the accretion of gas into 
the disc is likely to proceed through a Vi rial shock, and not 
along filaments dDekel & Birnboim 120061) . Our conclusions 
should therefore remain valid, so that massive galaxies should 



not provide winds powerfull enough to enrich the IGM by 
SNe feedback (but other way of enrich ment like gas strippin g 
from satellites can be considered like in ISchindler et ai]l2005l) . 
Dekel & Birnboiml d2006l) have computed the critical mass be- 
tween filamentary accretion and Virial shocks in the range 
10 11 -10 12 M depending on the metallicity of the infalling gas. 
Combining their results with ours, we can conclude that halos 
larger than a few 10 11 M Q cannot produce galactic winds, even 
within a realistic cosmological environment. At higher redshift, 
accretion rates are believed to be higher than the one we have 
considered here, so that galactic winds at earlier epochs are 
even more unlikely. A proper modelling of starburst (yet to be 
invented) might also provide an easier route for increasing the 
feedback efficiency of supernovae-driven outflows. 
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Appendix A: Numerical scheme for Supernovae 

We define the mass vanished in the star formation process as : 

{ Am s) SF = OT *(! + + rj w ) , (A.l) 

where the second term of the right part is equal to the total de- 
bris mass nid- We have also introduce the mass loading factor 
r]w, the parameter that determines the gas mass carried in de- 
bris. And m, is equal to 



1 + rjsN + W 



(A.2) 



6x is the volume of the cell where the star particle is spawned. 
Then if we introduce the fraction of debris particles : 



1 + t]sn + J] W 



(A3) 



we can simplify equation (IA.U by saying that a fraction fy of 
the gas consumed goes into debris particles, and a fraction 1 —fy 
goes into star particle. 

The maximum speed of a debris is given by 



I'd 



U-SN 



\f\ + rjwjj], 



(A.4) 



'.ViV 



us n is the typical velocity corresponding to the kinetic en- 
ergy released in one single supernova explosion (usn - 
3200 km.s- 1 ). 

Thus, the energy released to the gas by the debris is 



Ed - T]S1< 



Ms 



(A.5) 



Ms at and E$n are respectively the typical progenitor mass and 
energy of an exploding type II supernova (i. e. M SN = 10 M 
and E S n = 10 51 erg). The energy is independent from the mass 
loading factor. Then, we have to compute the moment conser- 
vation for constraining this free parameter. We want to repro- 
duce a Sedov blast wave, and so the velocity of the gas when 
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the debris are coupled with it is the Sedov speed of the shock 
propagation. 



1/2 



(A.6) 



/? is a multiplication factor nearly equal 1, we admit it is equal 
to 1, po is the density of the gas where the explosion takes place 
and Ax is the radius of the shock form the center of the explo- 
sion. We can also express Usedov in terms of u S n : 



U-Sedov — 



V2 



1 



(Sx\ 3 
\AxJ l+7] SN + r] W 



1/2 



USN 



(A.7) 



The Sedov blast wave is only valid in an homogeneous 
fluid, then it seems surprising to apply it in regions of star for- 
mation where the gas condenses in high overdensity regions of 
molecular clouds. That is the reality, but our simulations do not 
reproduce such high resolution features, and the assessment on 
the polytropic trend of the gas at high densities implies more 
pressurized regions in the disc. This behavior at high densi- 
ties transform high over-density regions of star formation into 
smoother regions at the scale of one cell of the numerical grid. 
That is the reason that allows us to use the Sedov blast wave 
model. We want the gas to have the moment of a Sedov blast 
wave, so the conservation of the moment gives us : 

m d u d - (m g + md)u Sed ov , (A.8) 

where m g = 4/3npoAx 3 is the gas mass where the supernovae 
explosion occur, and usedov — usedov/3 is the mean velocity of 
the gas carried along the shock. 

As we know all the quantities of this last equation, we can 
constraint our free parameter fd, and obtain 



1 4 /Ax 
1 + n\ — 

f d 3 \Sx 



V2 
15 



*(£) 



1/2 



= 1 



(A.9) 



If we want that our debris propagates to the nearest cells, let 
us assume that Ax 1.5, and we find the optimal value for 
fd — 0.5. Thus for a typical value of tjsn = 0.1 we find for the 
mass loading factor, our only free parameter, rjw - 1 ■ 
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Fig. 13. Mean velocity (up), velocity dispersion (middle), sound speed (bottom) of the gas as a function of the projected height 
for the Sd simulation (left) and the Ld simulation (right) at different times: t — 3 Gyr (solid line), t = 4.5 Gyr (dotted line), 
t — 6 Gyr (dashed line). 



